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Abstract 

We simulate the collision of atomic clusters with a weakly attractive surface using molecular 
dynamics in a regime between soft-landing and fragmentation, where the cluster undergoes large 
deformation but remains intact. As a function of incident kinetic energy, we find a transition 
from adhesion to reflection at low kinetic energies. We also identify a second adhesive regime 
at intermediate kinetic energies, where strong deformation of the cluster leads to an increase in 
contact area and adhesive energy. 



1 



The quantitative discourse on inelastic collisions between solid bodies stretches back at 
least as far as Newton's treatment in the Principia. Indeed, the inelasticity of collisions in 
flowing granular materials is frequently still described by Newton's coefficient of restitution, 
e, which is the ratio of the reflected to incident velocity (see Fig.[T]). Although the coefficient 
of restitution is often regarded as a material constant, generally it is found to depend on both 
the incident velocity and the adhesion between the solid objects^. Indeed, a full description 
of inelastic collisions is still a matter of on-going interest for the study of granular materials 2 . 
Nonetheless while inelastic collisions have been studied extensively for collisions of micro- 
or milliscale particles^^, much less is known about the collision of nanoscale objects such 
as atomic clusters. 

The collision of clusters with surfaces has been studied by many groups in the context of 
cluster deposition (for a recent review see Ref .-) . However, although the possibility of cluster 
reflection from hard surfaces has been considered (see the "phase" diagrams in Refs.- 1 ^), early 
studies focused on implantation^, fragmentation*^ or soft-landing^. In the so-called soft- 
landing regime, the incident kinetic energy is generally insufficient to overcome the adhesion 
between the surface and the cluster, resulting in collisions that always lead to adhesion. 
At high kinetic energies, one enters the fragmentation regime where the cluster fragments 
upon impact into atoms or several large components^, or the implantation regime where 
the cluster buries itself in the surface^ 2 -. Only recently has an intermediate regime been 
identified where clusters were observed to undergo a transition from adhesion to reflection^. 
This transition was exploited to assemble nanowires. 

Here we report on a detailed molecular dynamics study of the collisions of clusters, with 
kinetic energies that lie between the soft-landing and fragmentation regimes, on surfaces with 
a weak attraction to the cluster. The collisions span a range from weakly to strongly inelastic, 
and can lead to substantial deformation of the cluster. We have considered collisions of 
147, 309 and 561-atom Mackay icosahedra with a fee (lll)-terminated slab. The interaction 
between cluster atoms and between surface atoms was modeled using a single Lennard- Jones 
potential (in what follows a denotes the Lennard- Jones radius, e denotes the Lennard- Jones 
bond energy and m denotes the atomic mass). However the adhesion energy between the 
cluster and the surface was controlled by reducing the attractive (cr/r) 6 term between cluster 
atoms and surface atoms^ by a factor c (which was varied between 0.3 and 0.7). For instance 
when a 147-atom cluster was placed on the surface, it relaxed to form a contact angle of 
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120° at c = 0.7 and 100° at c = 0.5. In Re£2 a contact angle of 120° for Sb clusters on 
SiC>2 was estimated. Here we report results for c = 0.35 and only comment on trends for 
other values of c. We have also used a large cut-off distance of 6 a in the potential, which 
ensured that the adhesion energy between the cluster and the surface was well converged 
and enabled us to probe a collision regime where the kinetic energy and adhesion energy 
were delicately balanced. 

In each simulation, a central region of atoms in the slab follow Newtonian dynamics, an 
exterior 'shell' follow Langevin dynamics^ at a temperature T and a lower layer is fixed. 
In the Langevin region the friction parameter was varied linearly from at the Langevin- 
Newtonian interface to 2 at the Langevin exterior to minimize reflections of elastic waves 
from boundaries^. This shell of Langevin atoms regulates the temperature of the central 
Newtonian region and absorbs energy from the cluster impact allowing us to use a much 
smaller slab than would otherwise be possible. The size of the Newtonian and Langevin 
regions was chosen by studying a series of slabs of increasing depth and breadth until the 
energetics of test collisions were observed to have converged. The slab finally chosen was 11.7 
a by 11.3 a in plane and 10.3 a deep, comprising 5846 Langevin atoms plus 1344 Newtonian 
atoms. 

Figures [TH show the results from two collisions involving a 147-atom icosahedral cluster 
at an initial velocity v = 1.0 {ejm) 1 ! 2 with c = 0.35. Prior to collision, the cluster was 
equilibrated at a temperature of 0.13 e/ks for 10 4 time-steps and the surface was equilibrated 
at 0.2 e/kB (the cluster was equilibrated at a lower temperature to avoid evaporation). 
Figure [1] shows the velocity and deformation of a 147-atom cluster for a collision that results 
in adhesion. The potential energy of the cluster during this collision is shown in Figure [21 
Figure [3] shows the potential energy for an event that leads to a reflection. Note that in 
the collision that leads to adhesion, the cluster is edge-on to the surface. However in the 
collision that results in reflection, the cluster is orientated with a vertex closest to the surface. 
This difference in orientation leads to substantially different deformation of the cluster and 
dissipation of kinetic energy upon impact. 

Figure [2] shows both the change in internal energy of the cluster per atom AE C = E c — E^ 
(where Ef is the potential energy per atom of the isolated cluster) and the change in total 
potential energy per atom of the cluster AE pot = AE C + AE CS (where E cs is the cluster- 
surface interaction energy per cluster atom). At impact AE C rises sharply as the cluster 
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undergoes both reversible (elastic) and irreversible (plastic) deformation. Note that the 
loading and unloading phases of the collision are evident in both the change in the radius 
of gyration R g and the two-step change in the thermal kinetic energy E^ ermal . At the end 
of the unloading phase, we can see that while some of the internal energy of the cluster has 
been recovered (AE^ lastic ), most of the change in energy is irreversible ( AEp lastic ) . In the 
case of reflection (Fig. [3]) the cluster undergoes substantially less deformation than in the 
previous case, and the majority of the change in cluster internal energy is elastic. 

It is evident that the outcome of a collision depends strongly on orientation and thus 
cluster structure. However in what follows, we have averaged out this effect by conducting 
50 trials for each set of collision parameters, with cluster orientation varied randomly be- 
tween trials. We can then discuss the probability of bouncing at a given velocity and for a 
prescribed surface and cluster. 

Figure H] shows the probability of adhesion as a function of velocity. Strikingly, the 
probability of adhesion is a bimodal function of velocity at all cluster sizes with a second 
peak in adhesion at vq ~ 1.2 (e/m) 1 ^ 2 . At low velocities the cluster always adheres to 
the surface but at intermediate velocities the cluster starts bouncing. At higher velocities 
(0.5 < Vq < 3.0 (e/m) 1 / 2 ), we see reentrant adhesion followed by reflection in the limit 
where the velocity approaches the fragmentation regime. We find that this bimodality is 
also evident at c = 0.3 and 0.4, but disappears above c = 0.5, where cluster collisions almost 
always result in adhesion (at least up to fragmentation), and below c = 0.2, where only very 
slow collisions lead to adhesion. 

The average coefficient of restitution for each cluster size is shown as a function of initial 
velocity Vq in Fig. [5j Here we define the the coefficient of restitution as the ratio of the 
maximum velocity after the collision, Vf, to the maximum velocity prior to the collision, Vi 
(see Fig. [1]). Note that e is a relatively weak function of cluster size. At low velocities e 
is approximately constant but at high velocities, Vq > 0.7 (e/m) 1 ^ 2 , the coefficient exhibits 
a strong dependence on velocity: e ~ Vq ' 52 . This is a much stronger dependence on 
velocity than that given by Hertzian contact mechanics^ but is consistent with finite-element 
simulations of strongly plastic collisions with no adhesion 17 . We note that this dependence 
on Vq weakens in our simulations as c increases. 

Figure [5] also shows the relative change in radius of gyration, AR g /R g , of the 147-atom 
cluster at the moment of peak reflection velocity, Vf, as a function of vq. The fit shown is a 
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quadratic in v i.e. AR g /R g ~ Vq for v > 0.5 (e/m) 1 / 2 . Thus the relative spreading of the 
cluster during the impact on the surface, AR g /R g , is proportional to the incident kinetic 
energy. We note that this proportionality is consistent with strong plastic deformation, where 
the kinetic energy is dissipated largely at the cluster yield stress, Y, so that the plastic work 
~ YA(R 3 ) ~ YR 2 AR is proportional to the translational kinetic energy ~ VqR?. We have 
found that this relationship is relatively insensitive to both the values of c and the cluster 
sizes examined here. 

The amount of deformation affects the adhesion energy of the cluster by altering the 
contact area of the cluster with the surface. Fig. [6] shows the dependence of the adhesion 
energy Ei = —ET at the moment of peak reflected velocity Vf for the 147- atom cluster. For 
Vq < 0.5 (e/m) 1 / 2 , with little deformation of the cluster during the collision, Ej is constant. 
At intermediate to large velocities, as the cluster begins to deform on impact, this energy 
depends on velocity as Vq 50 . 

Also shown in Fig. [6] is the Weber number, We = E^/Ej, at the moment of peak reflected 
velocity. As the average value of We approaches and then exceeds 1, the number of clusters 
being reflected dramatically increases. In the low deformation regime (vq < 0.5 (e/m) 1 ^ 2 ), 
both the coefficient of restitution and the adhesion energy are approximately constant, so 
We ~ Vq. This is essentially equivalent to the liquid droplet models used in Ref^ to explain 
the transition from adhesion to reflection in their experiments. However, our estimate of 
the velocity at which strong plastic deformation occurs 5 , ~ 0.1 (Y/p) 1/2 -5-20 ms , 
is exceeded by the estimated cluster impact velocities of 100-200 ms" 1 in Ref.—. Thus it is 
probable that the collisions in Ref.— occur in the large deformation (plastic) regime. We 
note that the metallic clusters used in the experiments are likely to deform plastically at 
lower energies (relative to their binding energies) than the clusters described by the pair 
potential here^. 

At velocities Vq > 0.5 (e/m) 1 / 2 , both the coefficient of restitution and the adhesion energy 
become dependent on velocity as the cluster undergoes substantial plastic deformation. In 
particular, the reflected kinetic energy E^ goes as e 2 v^ ~ v® M for v > 0.7 (e/m) 1 ^ 2 and 
the adhesion energy E a goes as fo' 50 f° r v o > 0.5 (e/m) 1 ^ 2 (as shown in Fig. [6]). At the 
onset of this plastic deformation regime, E a grows faster than Ej , leading to a decrease 
in We and an increase in adhesion probability. Put simply, the increase in adhesion for 
v > 0.5 (e/m) 1 / 2 occurs because the clusters are highly deformed (they "pancake"), and 
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so there is an increase in contact area and (consequently) adhesion energy. However at 
velocities vq > 1.2 (e/m) 1 ^ 2 , we find that We increases as ~ v® A6 as Ej begins to dominate, 
leading to a corresponding decrease in adhesion probability. 

To summarise, we have identified a bimodality in the adhesion probability for clusters 
incident on a surface with poor adhesive properties. The re-entrant transition seen here is 
likely to be an example of a more general phenomenon which will occur when ever the Weber 
number at the transition between the low and strong deformation regimes is near one. Indeed 
as the Weber number depends on size, we expect that for a given cluster-substrate system 
there will be a range of cluster sizes where this re-entrant transition may be observed. For 
example, in metallic clusters, where the onset of strong deformation occurs at lower collision 
energies, we would expect the reentrant transition to be shifted to larger cluster sizes. 

We also find that the dependence of the coefficient of restitution on velocity is much 
stronger than that predicted by the classical theory of solid-solid collisions^ (leading to a 
strong dependence of the adhesion energy on velocity) although it appears to be consistent 
with finite-element simulations of strongly plastic collisions in larger particles^. This work 
suggests the possibility of experimental tests of the theory of collisions in nanoscale systems, 
in addition to the existence of new deposition regimes that could be exploited for device 
manufacture 1 - 3 . Further work will focus on the effects of cluster orientation and non-normal 
incidence on the cluster reflection. 
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FIG. 1: The evolution of the velocity and radius of gyration of the 147-atom cluster is shown over 
the duration of the collision illustrated in a series of snap-shots. The cluster has an initial velocity 
of 1.0 (e/m) 1 / 2 . In this case the cluster is captured by the surface and the velocity oscillates about 
zero after the collision. 
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FIG. 2: The change in cluster potential energy (both the cluster internal energy AE C and cluster 
total potential energy AE pot ) and the thermal kinetic energy (E^ erma i) are shown during the 
collision in Fig. (TJ 
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FIG. 3: The change in the cluster potential energy per atom (AE C ), the change in total potential 
energy per atom (AE pot ), the thermal kinetic energy per atom E^ ermal and the center of mass 
velocity are shown as a function of time in a collision that results in reflection. The cluster is a 
147-atom icosahedron and its initial velocity is 1.0 (e/m) 1//2 . 




FIG. 4: The probability of sticking as a function of incident velocity for each cluster size. At each 
velocity we conducted fifty trial simulations with randomly selected cluster orientations. Note the 
minimum and maximum in adhesion probability (at 1 and 2 respectively). 




FIG. 5: The coefficient of restitution as a function of initial velocity vo for c = 0.35 for all three 
cluster sizes. Also shown (on the right hand axis) is the relative deformation of the 147-atom 
cluster as a function of velocity which is well- fitted by a quadratic in vq. 
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FIG. 6: The adhesion energy E a at the peak reflected velocity and the corresponding Weber 
number as a function of vq for the 147-atom cluster with c = 0.35. The fit to the E a points is a 
linear function of wf]' 50 and is shown by the solid line. Note the maximum (at 1) and minimum 
(at 2) in the Weber number which correspond respectively to the minimum and maximum in the 
adhesion probability in Fig. [H 



